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1,0 INTRODUCTION 


A stochastic simulation .ncdel of an Optical Contrast Seeker (OCS) was developed and 
incorpora‘cd into the existent Monte Carlo point target, laser guided, termina) homing 
6-DOF simulation program described in reference 2, Two OCS simulation subroutines 
were developed: (1) a subroutine consisting of high and low frequency poles, and (2) a 
simplificetion of the first consisting of only low frequency poles. The non-essential high 
frequency poles of the OCS model were eliminated to reduce computer run time. This 
simplified subroutine is, In general, sufficient; however, either subroutine may be used by 


setting the appropriate flag. 


The 6-DOF si.7ulation program was modified to inchide a two-dimensional target and 
helicopter vibrations, The target dimensions are required for study of the OCS breaklock 
and blind-range phenomena, The helicopter vibrations were added to describe the missile 


launch transients. 


The updated 6-DOF simulation program may continue to be run as either a stochastic or 
deterministic program. In the stochastic mode, the program executes a specified number 
of runs*, computes miss distance coordinates from each run, and then determines the CEP 
and other statistical parameic:s trom the set of miss distances. Each run of the run set is 
made based on both initialization error conditions (mass unbalance, etc.) and time varying 
error conditions (wind, etc.) that are randomly generated from input error probability 
distributions, In the absence of statistical input data, the operation of the 6-DOF program 


reverts to that of the deterministic version of the program. 


Because of the addition of the two OCS subroutines and associated 6-DOF program modifications, 
the computer program listing contained in this document supersedes the listings found in 
references 1 and 2. However, the deterministic program model description and input/output 
formats found in reference 1 and the stochastic program description and input/output formats 


found in reference 2 are still valid. 


*A run is defined as the numerical integration cf a missile trajectory from launch to 
target plane intercept. 
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2.0 FUNC TIONAL DESCRIPTION 


Figures 2,1 and 2, 2 are the block diagram representations of the OCS simulation model, 
Model S2, Figure 2.1 represents the pitch channel and Figure 2.2 represents the yaw 
channel. The transfer functions given in these block diagrams were transformed by use of 

the M-method into state variable format for solution by numerical integration. A aew seeker 
subroutine (S2) was developed to control the integratior of these state variables and to perform 
other cilculations pertinent to the OCS, 


Because of the modular concept of the 6~DOF simulation program, interfacing the new sub- 
routine required minimal effort, particularly since the input and output variables of this 
seeker subroutine were identical to those of the laser seeker subroutine, Sl. Tho one 
exception to this was the missile-body-to-seeker-gimbal coordinate transformation matrix, 
since the OCS sceker mounting has been rotated 90 degrees from that or ‘he laser seeker. 
This necessitated a new matrix for transformation of the line-of-sight vector ;LOS) from 
missile body coordinates (the basic program coordinate system) to the seeker gimbal 
coordinate sy stem. 


Figures 2.3 and 2.4 are the block diagrams of the simplified model, Model S88, containing 

only the low frequency components with the high frequency tiocks replaced by their appropriate 
gains. This simplification, in effect, neglects the small amplitade high frecuency oscillations 

of the system which are superimposed on the more significant, lower frequency dynamics. 

The simplification e ived as a consequence of the very long computation time required when 
running Model S2. ‘Tne high frequencies of 314 and 1000 rad/se« of Model S2 require a numerical 
integration step size of .5 millisec, wh‘le Model S3 is integ:-ated accurately at a step size of 
12.5 millisec. This translates directly into a factor of 25 difference ir. ~un time, anc: there 

is no significant loss of accuracy when using Mode! $3. In comparison runs that were made, 


less than three percent differences in miss distance and CEP were observed. 
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: 2.1 Transformation Matrix i 


The seeker gimbal coordinate system Xe YG » Z g) is given in Figure 2.5 with respect to 
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the missile body covrdinate system > Yp> Z,)- The order of rotation is (1) a rotation i 

about Y, through the outer gimbal angle @., then (2) a rotation about Z G through the inner i 

gimbal angle . 
The missile body-to-seeker gimbal coordinate system transformation matrix for this rotation 
t sequence is: 


me 9g vos tq) sin YG (-cos Yo sin 8g 
[M], = |(-cos 6, sin Yq) 98 (sin Yq sin 8q) Eq. 2.1 | 


sin 8G 0 cos 8c 


(1) bawst Rotation About 7, B 


(3) Missile Body/Gimbal Coordinate System 
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FIGURE 2.5. Coordinate Rotation 


SLE TET ET HT 


2 pe REP ERMNT P ERTS EROS 7 CHOTA ERY SPREE ES OF TET 
AE PERE SLE EE NE FE FT PE I EI ET eT I ETS RET IIT I SST a - 


2,2 OCS Support Mocels 


In addition to incorporation of the pitch and yaw axis OCS models into the 6-DOF program, 
other critical parameter models unique to the OCS were added. These parameter models 


are listed below and described in the following sections: 


@ Gimbal friction 

e Target size 

e Seeker breaklock 

e Seeker blind range 

r) Transport lag 

e Helicopter induced launch transients 
2.2.1 Gimbal Friction Model 


Gimbal friction couples the missile angular rates to the OCS platform and may cause system 
degradation. Gimbal bearing pre-load and the gimbal torque motor are the main contributors 
of this friction which is primarily sticticn (static) and coulomb. Coulomb friction is defined 
as a constant frictional drag which opposes motion but has a magnitude that is independent 

of velocity. A slight disunction must be made between stiction and coulomb friction because, 
in general, the force required to initiate motion (overcome stiction) is somewhat greater 
than the coulomb friction. When the stiction level is identical to the coulomb friction and the 
system starts at rest, any applied force to the gimbal less than this value must be identically 
opposed so that no motion is initiated. Thus, the idealized coulomb friction model, 

Ty = To sen(6,)s can create a physically impossible situation where the friction model 
supplies energy to the system. As mentioned previously, the missile rates (r. -ually, 
accelerations or torques) are directly coupled to the platform through the effects of friction, 


although the coupling is limited in magnitude by the friction level. 
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Figure 2.6 is a block diagram revresentation of the pitch gimbal friction model. Tlie opposing 
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friction, Ty is a function of the coulomb friction, TG} the net torque, Ty and the relative 
gimbal rate, 85 (Stiction is assumed identical to the coulomb magnitude.) The missile angular 
acceleration, Q, is coupled to the platform through the function fg which is magnitude limited. 


The complete mathematical description is: 


For| édl< A, where 4 is a computational dead zone and the system is considered at rest, 


Ty = Ty Sat(T.,) Eq. 2.2 

ToT. Sen(Q) 

Qa, =Q Sat = » where J is gimbal inertia. Eq. 2.3 
For ig 4 

T= T, Sen(,) Eq. 2.4 

Qn = 0 whore Q,d i, (Tg Ty Q, 6a Eq. 2.5 
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FIGURE 2.6. GIMBAL FRICTION MODEL 
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Figure 2.7 is a time history of the torque motor output. For this illustration, 8c is 
initially zero ard 8c is positive. Te initiate motion, the torque tnotor had to exceed the 
coulomb friction level and the wire compliauce torque. Figure 2.8 shows the missile 
acceleration, Q, and the gimbal acceleration, Ws Observe that the gimbal angular 
acceleration tracks the missile angular acceleration until the torque motor exceeded 


the breakaway torque level. 
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Figure 2.7. Pitch Channel Torque Motor Output 
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Figure 2.8. Missile Pitch and OCS Outer Gimbal 
Inertial Angular Acceleration 
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2.2.2 Target Model 


A two dimensional target model was incorporated into ths program for the purpose of 
computing seeker breaklock (Section 2. 2.3) and blind range (Section 2.2.4). The 


dimensionality of the target does not affect any other missile or secker parameter. 


The target model is defined by the shaded area in Figure 2.9 and due to the dimensionality 

is restricted tu be normal to the line-of-sight of the missile, The outer rectangle represents 
the seeker field-of-view (FOV) and the intersection of the dashed lines the instantanecus 
seeker aimpoint. The top and bottom of the target are always parallel to the raster lines 

of the vidicon screen of the seeker. Thus, for this elementary model, if the seeker roils, 


the target rolls. The line-of-sight vector intersects the target ai its geometric cen‘er, 


FOV 
(So pe OUTLINE 


FIGURE 2.9. Target Model 
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2.2.3 Seeker Breaklock 


Seeker breaklock (loss of target image) is assumed* to occur when sceker aimpvuint shifts 
more than 50 percent of target height or width during one sample period, 7. Seeker aimpoint 
is defined as the projection of seeker boresight onto the target nlane. Fignre 2.10 illustrates 


the geometry of the breaklock parameter, aimpoint shift. 


S.. is the total amount of seeker aimpoint shift over the sample period 7. S,, and Sw are the 


H 
components of 8, that parallel target height Ty and width Ty respectively. If Sy becomes i 


greater than .5 T,, or if Sw becomes greater than .5 T,,,, it is assumed that breaklock has 


H 
occurred. When breaklock does occur, the simulation run is terminated because breaklock 


causes loss of the missile as far as homing in on the target is concerned. 


Seeker aimpoint 
at time =t 


Seeker aimpoint 
at time =t-<7 


Target 


FIGURE 2.10, Seeker Aimpoint Shift 


*This assumption is valid for a target-to-range ratio as low as .001. 
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In the simulation program, breakiock is determined from the angles subtended by the target 


sides/slant range and ain-point shift/slant range. Thus, for pitch, the angle* 


ee | 
Ore = tan [(T,,/2)/R] Eq. 2.6 
and the angle 
be 7.6 -B ) Eq. 2.7 
Z Bz Zit 


are compared at the end of each sample period, T. If 4, becomes greater than 


3) breaklock conditions are met. 


TH’ 
The variables in the above equation are defined as: 


R ~ slant range 


G14,  ~ angle subten‘led by slant range and one half the target height 


By -~ angle subtended by pitch plane component of LOS S)) and seeker 
ter boresight, one sample period back in time 


Br, = angle subtended by pitch plane component of LC'S S) and seeker 
: boresight, at the current time 


Ay — angular shift of aimpoint over one sample peyiod 


When breaklock is encountered, the message be!ow is printed out aud the run is terminated, 


The time (seconds), range (feet from targ-t), and tt channel in which breaklock occurs is 
output. 
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An 2xample of seeker aimpoint shift is given in Figure 2.11, This figure contains a time 
history of seeker aimpoint shift for a 1 kilometer deterministic trajectory. No errors (such 


as gimbal mass unbalance, launch transients, rate gyro drift, etc.) were present when this 
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run Was made. The sample period was 16.7 millisec. The breaklock point (50 percent of 


*Comparable calculations are made in the yaw plane to determine if breaklock occurs there. 


le 
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target) currently in use in the simulation program is illustrated in the figure. If any of the 
spike’; on the curve had reached the 50 percent point, the trajectory would have been 


terminated. As noted for this particular ex: mple, the maximum boresight shift was about 
ten percent. 
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FIGURE 2.11. Seeker Breaklock Time History 
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2.2.4 Seeker Blind Range 


' Secker blind-range is defined as the distance from the target when the target image occupies 
70-80 percent of the tracker vidicon field-of-view (FOV). The actual blind-range limit is 
adjustable and is currently preset for 70 percent FOV. 


To determine if a blind-rai.,e condition occurs, gating functions are established to check 

the target edges with respect to the 70 percent FOV lines as shown in Figure 2.12. However, 
with this implementrtion, biind-range is a function of the LUS angle, as well as target growth. 
Figure 2.13 is a tine histcry of target edge in one axis with respect to the vidicon FOV. The 
general shape of the curve is ue to target growth while the perturbetions are caused by the 
LOS angle variations. The piogram monitors all four edges at eaci. integration step to 
determine if any one edge “4 reached the blind-range limit. Noriinally, this 'imit is set at 
70 percent FOV, butisa ae tram variable that may be input at any desired value. When the 
blind-range limit has been rzached, the seeker rate gyro's output signals to the autopizot are 


| 
[ 


held at their present value. Thus, the missile no longer responds to commands gencrated by 
the seeker, but instead, flies into the target with the av pilot signals set at the blind-range 
value. In Figure 2.18, blird-range occurred at 3.61 x “nds. When blind-range occurs, the 
‘ullowing message is printed out: . 


ek OR OK OR Ok OR Ok Ok RR RK RR Ok Rk OK KK OR OK OR OK OK KOR OK OR OR OR RK OR OK OK OR OK OK OK 


CCS BLIND RANGE SIGNAL HOLD AT TIME = 3.01 RANGE = 225. 00 
KOR OR OR ROR KOR KOR EK OK OR OR OK OK KR KOR OK OK OKO OR RR OR OR OK KOK OK OK Ok OR RO OK OK OK 


Time is in seco id» of flight and range is in feet irom target. 


VID.CON ' BLIND-RANGE 
¥OV o ‘ GATING FUNCTIONS 
' 
BEN eee - : 
i] 


BORE 
SIGHT 


FIGURE 2.12. Blind-Range Due to LOS and Target Growth 
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2.2.5 Transport Lag Implementation 


The exponential, e SF in the second block of each of the block diagrams of Figures 2.1, 


ihethgsln dc Salam GA iy Gat Fy tt Ba se ok ne 


2.2, 2.3, and 2.4 represents u transport lag in the camera of the OCS, The actual phenomena 
that this term models is unknown at the present; however, it closely approximates the effects 


required to match OCS subsystem test data. Implementation of this lag was accomplished 
as fullows: 


Jeena soci Bb 


tae hes DRS Re Moe OB TSB att on 


® The last six time points of the LOS angle (BZ and BY) are stored in 
a storage array. 


e A table look up function interpolates linearly within the stored values to 
return the values of BZ and BY at the laz time (TIME - TLAG), 


e Tre interpolated values of the LOS angles are then used by the OCS model 
as the target position fur guidance calculations, 
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2.2.6 Launch Transients 


: Missie pitch, yaw, 2nd roll launch transients caused by helicopter vibration and launcher 
rail/missile shoe interaction were modeled because tnes2 iransients could severely affect 
the condition of seeker breaklock. Typical launch trausient data from which the models were 
developed are given in reference 3. Examination of this data reveals that pitch and yaw 
characteristics are similar. Thus, pitch and yaw perturbations are simulated from similar 
transient models, Roll transients are presented as clots of roll angle versus time. Roll 

: transients appear to exhibit diffzrent characteristics from those of pitch and yaw, thus 

roll is modeled separately, An cxample of the plots as presented in reference 3 is given 


in Figure 2.14. 
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FIGURE 2.14. Typical Launch Transients 
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In addition to the data of reference 3, telemetry data relating to pitch and yaw, in the form 
of undocumented strip charts for pitch and yaw rajes and a power spevtra) density (PSP) 
plot of rate-gyro output in pitch, were available to aid in model development. These 
telemetry data were taken from a captive flight test of the anti-tank missile mounted on a 
helicopter launcher. A reproduction of the PSD is shown in Figure 2.15. 


2.2.6.1 Pitch, Yaw, and Roll Models 
2.2.6.1.1 Pitch Model 


In modeling the pitch launch transients, the helicopter vibration was assumed to be a pitching 
moment applied to the launch rail and coupled to the missile through the rail shoes, Since 
the power spectral density (PSD) of the pitch rate gyro output, Figure 2.15, has distinct 
frequencies; the pitching moment was modeled as a harmonic forcing function and is applied 


until the rear shoe exits the launch rail (tip-off point), The forcing function is defined as: 


At yn 
Fit) = Ase ©. 5 A, sinyt+@,) Eq, 2.8 
m i i i 
irl 
where 
wo frequency (217 f.) 
a - phase angles initialized randomly from a 0 to 27 
uniform distribution 
A, - peak amplitudes 
A, - time constant which spreads the PSD about 
the discrete frequencies 
A - scale factor 
m 


The relationships between the peak amplitudes, Ay were determined by comparing the peak 
amplitudes of the PSD at the desired frequencies. However, these values are the pitch rate 
densities in (rad/sec)”/hz and need to be related to the equivalent autocorrelation values. 


Since the autocorrelation function, R(0) is related to the PSD by 


R(0) = f PSD (f) df 
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It can be shown that with an ideal bandpass filter of bandwidth, BW, we car relate the 
autocorrelation values, A,, to the PSD values, By by 


Using the first four (n = 4) predominant frequencies of the PSD given in Figure 2.15 


f, = {11, 22, 33, 44} 
and their peak values 


B, = {.085, 0.10, 0,56, 1.0} 
the autocorrelation values are calculated with BW = 500 hz as 
A, = {6.5, 6.7, 16.7, 22.4} 


Usirg the time constant, A, =-1, to spread the spectral densitics about each frequency, the 
* 
scale factor, An was determined by calculating the autocorrelation function, Figure 2.16, 


ef the simulation rate gyro output signal and comparing it to the telemetry PSD, Figure 2.15. 
A reasonable match was obtained with all the telemetry data with 


A, itch) = 10 FT-LBS 
and 


A, = {1, 4, 12, 26}. 


Example time histories of F(t), the pitch rate, and the rate gyro output are shownin Figuce 2.17, 


2.18, and 2.19, respectively. The missile was constrained to the launch rail for these runs. 


An example of a pitch rate time history of the missile (not constrained to the rail) for the first 
one-half second of flight is given in Figure 2.20. Front shoe exit time is .086 seconds and 


rear shoe exit time (tip-off time) is .112 seconds. 
2.2.6.1.2. Yaw Model 


The pitching moment forcing function and coefficients were assumed applicable for yaw. The 
peak amplitude, A of yaw moment was determined in similar manner as pitch. Amn (yaw) 
was adjusted to give a reasonable match between the peak-to-peak yaw rates of the simulation 


program and the peak-to-peak rates given in reference 3. An example of a yaw rate time 


* The autocorrelation function was calculated with the Time Series Analysis Program 


(TSAP), Reference 4. 
20 
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Figure 2.16. Autocorrelation Function of Rate 
Gyro Output - Simulation Model 
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FIGURE 2.17. Moment Forcing Function Time History 
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FIGURE 2.18. Pitch Rate Time History 
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FIGURE 2.19. Output to Autopilot Time History 
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history is given in Figure 2.21. The best value obtained for the yaw rate amplitude 


coefficient was: 
An (yaw) = 4 FT-LBS. 


2.2.6.1.3. Pitch and Yaw Rate Initialization 


As long as the missile is on the launch rail and the pitch and yaw moments are being applied, 
the pitch and yaw rates are directly proportional to the integral of F(t). Integrating F(t) 
results in the following equation, H(t). 


At Ay A 
ad e i 1. A 
H(t) = Awe r >, 2 rd [A, sin (at +) - a, cos t,t +,)) Eq. 2.9 
i=l e i 


The proportionality constants for pitch and yaw are the moments of inertia about the 
appropriate rotational axes. Thus, 


we = 2 Eq. 2.10 
We = ie: | Eq. 2.11 
Zz 
where 
H(t) Q 7 refers to the pitch rate equation 


refers to the yaw rate equation 


= 
Gs 
w 
' 


ly - moment of inertia about Y axis 
I, - moment of inertia about Z axis 


Equations 2%. 10 and 2.11 are solved at time = 0 to determine the correct initial values of 

pitch and yaw rate so that the time functions of pitch and yaw are in phase with the moment 

time function F(t). (This is equivalent te solving for the constant of integration of a differential 
equation.) Pitch and yaw rates are not determined by Eq. 2.10 and 2.11 excent for initialization. 
Instead, the forcing function F(t} is numerically integrated along with all other differential 
equations in the program. 
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FIGURE 2.20. Simulation Model Pitch Launch Transient 
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FIGURE 2.21. Simulation Model Yaw Launch Transient 
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2.2.6.1.4 Roll Model 
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Roll transients were modeled based solely on the data given in reference 3. No additional : 
data on roll was available. Ati roll data plotted in reference 3 is given in roll angle; however, 
as with pitch and yaw, the roll transient was modeled as a roll moment. A roll moment model 
was developed that produced roll angles that matched typical roll angle characteristics as 


given in the reference. * 


b Examination of the slopes of the plots of reference 3 reveal that for all practical purposes, 
roll rate is zero until the front shoe exits the launch rail, Roll rate then begins to build up 
until a maximum roll angle is reached between .3 and .4 seconds. Once off the rail, roll 
angle varies about a mean value established when the autopilot gyros were uncaged. The 


oscillation about this mean is caused by the roll stabilization system. 


The roll transient model developed assumes a zero roll rate and a zero roll acceleration (zero) 


moment applicd) until front shoe exits the launch rail, as shown in Figure 2.22. Following 


SRN = SNE SIE PTY PP AES Sno 


Yront shoe exit, a roll moment (F yx) is computed based on the difference between front and 
rear shoe exit times (Qt) and the desired rear shoe exit roll rate (tip-off roll rate, Woo)? 


or: 


Fux = Woro . 1,)/At Eq. 2.12 


where x is the moment of inertia about the longitudinal axis and wipro is specified by the 
user. The exit time difference (At) is .026 sec for this missile simulation, 


The roll moment F yx) is applied to the missile from the time of front shoe exit to the time 


of rear shoe exit. Integration of the roll acceleration due to F results in the tip-off roll 


MX 
rate Won’ Immediately following rear shoe exit, all transient models are zerocd out and 
the simulation reverts to 6~DOF guided flight. An example roll angle time history resulting 


from this model is shown in Figure 2. 23. 


Sebati A ce ga Ra 


* Modeling of roll transients is critical because of the interaction between seeker break- 
lock and roll acceleration. Seeker breaklock is influenced by roll acceleration through 
the roll coupling term of the seeker rate gyro output axis. This output axis is directed 
along the roll axis of the missile. 
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FIGURE 2,22, Launcher/Missile Configuration 
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FIGURE 2.23. Simulation Model Roll Launch Transient 


26 


: Shad male yeti ra 
See ty tas evade tina tbe RE a tee etn ne 
a pe ce cae we ok ee tenes ait ee a a ll A me eR 


eer er rene a Sa a i i 


F 
? 
E 
C 


ed 


E 
[ 
k 


2.3 Input Variable Description 
2.3.1 OCS Input Variable Description 


The following list identifies all variables of the OCS seeker subroutines (S2 and S3) that can 
be input by 3-cards.* Variable names beginning with K are the gains shown in the block 
diagrams of Figures 2.1 and 2.2. Variable names beginning with W are the frequency 
components given in Figures 2.1 and 2.2. The frequency variable names themselves do 
not appear in the figures. Instead, actual frequency values are shown. However, these 
variables may still be input by 3-cards. The variable name and its position in the block 


diagram can be correlated by the frequency values tabulated below in the second column. 


SYMBOL 
FORTRAN USED IN c 
SYMBOL TEXT INDEX 


DEFINITION 


* PITCH CHANNEL GAINS 


YAW CHANNEL GAINS 


Input description continued on next page. 


. See reference 1 for definition and use of 3-cards. 
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SYMBOL 
FORTRAN USED IN ; 
SYMBOL TEXT DEFINITION 


WTQl 
: WTQ2 
WGQl 
WGQS 


WGQ4 : PITCH CHANNEL FREQUENCIES 
WGO5 
V GQ6 
WRQ2 
WRQ4 


SOEUR TE, PLE FY FORTE TPE NIE TE ETS 


YAW CHANNEL FREQUENCIES 


Rate command limit in pitch und yaw 
Torque command lirnit in pitch 

Torque command limit in yaw 

Moment of inertia of inner gimbal 
Moment inertia of outer gimbal 

Rate gyro gain to autopilot, pitch and yaw 
Inner gimbal friction coefficient (in-oz) 
Outer gimbal friction coefficient (in-oz) 


Blind range decimal percent field of view 


Target height (ff) 


Target width (ft) 


cee ee Tahoe 


Seeker sample period (sec) 


OCS transport lag (sec) 


i 
| 


2.3.2 Launch Transient Input Variable Descrip’ion 


A new subroutine, LTRAN, was added that contains the pitch and yaw moment forcing 
function, F(t), and the pitch and yaw initialization function, H(t), This subroutine is called 
by subroutine A3i for initialization of the rates and by subroutine A2 to compute the time 
varying rates. 


The roll transient model is also initialized in subroutine A3I. Subroutine A2 contains the 
logic for integration of the roll acceleration computed in A3I. 


The following list identifies all variables of the launch transient models that are input by 


3-cards. 


SYMBOL 
USED IN c 
TEXT INDEX DEFINITION 


Tip-off roll rate (deg/sec) 


Peak amplitude of pitch moment forcing function 
(ft/lbs) 


Peak amplitude of yaw moment forcing function (ft/lbs) 


Launch transient vibration flag (pitch and yaw only) 


0 - no vibration 
1 - run With vibration 
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3.0 RANDOM ERROR SOURCES 


Initial condition random error sources specified as probability distributions unique to the OCS or 
impacting the operation of the OCS are listed Selow and described in the following sections. 


J, Seeker Platform Mass Unbalance 
e Cuter gimbal 
e Inner gimbal 


2. Seeker Rate Gyro Errors 


e Drift 
® Mass Unbalance 
© _ Output axis/missile roll coupling 


3. Launch Transient-Rate Distributions 


e Pitch and yaw rate 
e Roll rate 
3.1 Error Source Distribution 


3.1.1 Seeker Platform Mass Unbalance 


Error randomization of seeker platform mass unbalance was added to the simulation program 
for both inner and outer gimbals of the seeker head. Missile acceleration normal to the 
gimbal plane acts on this mass unbalance to create a torque that attempts to rotate the 


seeker head, thus generating an error signal. The seeker torque motor must then compensate 
for this error. 


Figure 3,24 illustrates the mass unbalance geometry of the outer gimbal ring. The centroid 
of mass unbalance is assumed to have an equal likelihood of lying at any point in the gimbal 
plane, while the mass unbalance magnitude Kio) distribution must be specified by the user. 
Since mass unbalance magnitude includes the moment arm as well as the normalized force 
acting on the moment arm, the actual radius on which the mass unbalance lies is unimportant, 
except as to whether the resulting torque causes a clockwise or counter-clockwise rotation 
about the gimbal axis. Thus, following the selection of the mass unbalarce magnitude from 

a specified distribution, the sign of the mass unbalance is selected from the uniform 


distribution shown in Figure 3,25. The torque acting about the outer gimbal axis is then: 
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FIGURE 3. 24, Qeeker Outer Gimbal Mass Unbalance 
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FIGURE 3.25. Probability Distribution of Mass Unbalance Sign 
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Typo ~ Axo’ Kuno * SGN(RN) Bq 
{ 
The total acceleration, 


Ar = [Ay p? Ayp Avy) Eq. 2.14 


acting on the missile is resolved through the outer gimbal angle (@q) to get the acceleration 
component (Aygo) normal to the gimbal plane. Thus, from Figure 3. 26, 


AxcGo = Ayn cos 8G - Arn sin 8c" Eq. 2.15 


Since the outer gimbal rotates about the y-axis ,) of the missile, only the acceleration 
components along missile x-axis (X,) and z-axis (25) act on the gimbal plane, 


FIGURE 3.26. Gimbal Rotation About Missile Y-Axis 
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The inner gimbal mass unbalance centroid is somewhat more complex to locate because the 


inner gimbal mass is more complex in shape as shown in Figure 3. 27. 


FIGURE 3.27, Inner Gimbal Mass 


As with the outer gimbal, the distribution of the mass unbalance magnitude Kopp is user 
specified, while the plane in which the centroid lies is assumed to have an equal likelihood of 
being oriented at any angle ( X ), as shown in Figure 3.28, The angle xX is selected from the 
uniform distribution shown in Figure 3,29. The resulting torque acting about the inner gimbal. 
axis is then, 


Tupr ~ Supr 4yg 08 X ~Axg sit X) mine 
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FIGURE 3. 29, Probability Distribution of Mass 
Unbalance Radial Position 
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The acceleration components (Ay G and Ay) in the inner gimbal coordinate system are 
determined by transforming the total acceleration A.) through the gimbal angles 
(8, and Yq) thus 


AxG AxR 
al = Us: Aye Eq. 2.17 
ATG Arn 


The transformation matrix [M] G is derived in Section 2.1. 


Mass unbalance magnitude distributions are normalized with respect to acceleration, Thus, 
mass unbalance has the units of IN-OZ per ¢g of acceleration. All acceleration components 


seen in the preceeding equations are also normalized. 
3.1.2 Seeker Rate Gyro Errors 


The seeker rate gyros can produce errors that perturb and bias the tracker rate loops and 
degrade the performance of the OCS seeker. Therefore, the primary errors sources of these 
rate gyros were included in the simulation program for study. Gyro errors (or more correctly, 
gyro error torques) can arise from a variety of different sources and are usually expressed as 


equivalent gyro drift rates. 


Constant gyro drift rates result from uncompensated bias torques, Ls and their magnitude 

is usually a measure of the gyro quality. Acceleration sensitivity is primarily a function of 
the mass unbalance about the gyro output axis, i.e., the center of mass not coincident with the 
output axis. Linear accelerations, A, normal to the output axis produce gyro drift rates 


proportional to the acceleration. 


The last gyro error source to be considered is a characteristic due to the inertia of the gyro 


float assembly. Angular accelerations about the output axis, » cause gyro pick-off angle 


Won 
errors, which in turn torque the gyro through the electronic caging loop resulting in a gyro 
drift rate. Thus, the total gyro error rate is 


Jon 


A+ Won Eq. 2.19 


re 
w= a + 


mI 
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where H is the angular momentumof the rotor, P is the pendulosity, and J OA is the moment 
of inertia of the floated assembly about its output axis. 


The OCS platform has two stabilizing rate gyros mounted to sense inertial rates in the pitch 
and yaw axes, These rate gyros are so mounted that their output axes are aligned with the 


platform roll axis. Therefore, from Eq. 2.19, the error rate for the pitch gyro is 


Weg = Kpo *Kpo A * Kono YRoLL Eq. 2.20 


Wey = Kay +KpyA + Koar Monn * Eq. 2.31 


The coefficients in Eqs. 2.20 and 2,21 are user specifiec by their distribution functions. 
Since the gyro pendulosity, Kp has an equal likelihood of occurring anywhere about the 
gyro output axis, its location is picked from a uniform distribution prior ‘ce each run, 
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fe 3.2 Probability Distribution Input Description 
: 3.2.1 OCS Monte Carlo Input Variables 
The variables associated with the Monte Carlo seeker models are given below. The mean 
F values of these variables are input by 3-cards and the probability distributions are input by 
: * : 
t 8-cards, 
Program Cc 
t Variable Index 
Name of of MCARLO Flag* 
c Ecror Error . 
b Source Source Definition 
KUO 611 S21 Ovtex gimbal mass unbalance 
(in-oz/g) 
KUI 612 S21 Imer gimbal mass unbalance 
(in-02z/g) 
KBO 613 S21 Outer gimbal drift rate 
(deg/sec) 
KB! 614 S21 . Inner gimbal drift rate 
. (deg/sec) 
KPO 615 S2I Outer gimbal pendulosity 
coefficient (deg/sec/g) 
KFI 616 S21 Inner gimbal adrift rate 
(deg/sec/g) 
KOAO 617 S82] Outer gimbal output axis/ 
roll coupling coefficient 
: (sec) 
KOAI 618 S821 Inner gimbal output axis/ 


roll coupling coefficient 
(sec) 


*MCARLO is flagged by the C-Index of this variable in the calling module, 


When MCAR:_O is flagged by this C-Index, a random number will be returned from 
MCARLO for the error source in the first column, 


ay 


See reference 2 for definition of 8-cards. 
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3.2.2 Launch Transient Monte Carlo Input Variables 


The variables associated with the Monte Carlo launch transient models are given below, An 
8-card is used to select any one of these models (roll, pitch or yaw) as a Monte Carlo variable. 
Roll is the only one of the three that requires specification of a probability distribution on the 
8-card, The pitch and yaw models do require 8-cards; however, the probability distribution 


input fields are left blank because pitch and yaw are randomized indirectly, as explained in 
Section 2, 2.6. 


A mean value of roll rate (WPTO) is input by 3-card. Mean values of pitch and yaw rate are not 
input, because the mean and distribution of these two variables are determined from solution 

of the forcing function, F(t), Eq. 2.8. However, the peak amplitude of pitch (AMP2) and yaw 
(AMP1) moment (due to helicopter vibration) must be input by 3-card. In addition, the flag, VIB, 
defined in Section 2, 3,2 must be input equal to 1, 


Program 
Module 
Calling 
MCARLO 


Ail, A2 Mean tip-off roll rate 
(deg/sec) 


A3I, A2 Peak amplitude of pitching 


moment forcing functlon 
(ft/Ibs) 


ASI, A2 “| Peak amplitude of yawing 
moment forcing function 
(ft/lbs) 


*MCARLO is flagged by the C-Index of this variable in the calling module, 


When MCARLO is flagged by this C-Index, a random number will be returned from 
MCARLO for the error source in the first column. 
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3.2.3 Pitch and Yaw Randomization Independent of Launch Transient Model 


Pitch and vaw tip-off rates may be randomized from an input probability distribution by 
inputting the C-indices of pitch and yaw rate on an 8-card. This capability was added as an 
option to directly randomize as opposed to indirectly randomizing pitch and yaw rates as 
mentioned above. Use of this option will generate instantaneous changes in pitch and yaw 
rate at time of rear shoe rail exit. This option was added primarily to allow randomization 


of pitch and yaw rates for launch from a tower cr ground vehicle in which there are no 


launcher vibrations. However, this option can be exercised simultaneously with the vibration 
model above. Rol! rate randomization described above applies equally to helicopter or 


ground launchers. 


Variable Program 
Name of Module MCARLO Flaz* 
Source MCARLO Name. Index Definition 


Se a: ak heat a i Ta al la ale et a re 


Pitch rate (deg/sec) 
Yaw rate (deg/sec) 


*MCARLO is flagged by the C-Index of this variable in the calling module, 


When MCARLO is flagged by this C-Index, a random number will be returned from 
MCARLO for the error source in the first column. 
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4.0 COMPUTER PROGRAM DESCRIPTIONS 
4.1 New Subroutines 
; 


The basic structure of the 6-DOF Monte Carlo program remains unchanged, Modifications 
i and minor alterations were made to incorporate the OSC subroutines and all related models. 
Four new subroutines were created and minor changes in other subroutines were made to 
interface the uew subroutines with existing program structure. The new subroutines added 


are: 


1. LTRAN - pitch and yaw helicopter vibration subroutine 
containing Eq. 2.8, the moment forcing function, 
and Eq. 2.9, the rate initialization function. 


2. $21 - OCS initialization subroutine. This subroutine 
initializes values for both ~CS models, S2 and S3. 
(S2I contains an entry point, S3I, that initializes 
variables pertinent to S3.) 


3. 82 - The derived OCS model subroutine, which contains 


both high and low frequency components. 


4, S3 - The simplified OCS model subroutine, consisting, of 


only the low frequency components. 


Input variables for the OCS model and all related models are defined following the sections 


that describe those particular models. 


Seeker subroutine selection is made by use of the 2-cards? Care must be exercised when 
changing 2-cards, because, as explained in reference 1, the order in which 2-cards are 


input determines tne order in which all missile and environment subroutines are called. 
4.2 Monte Carlo Runs With Breaklock 


Flights that breaklock are terminated when breaklock occurs. This, of course, impacts 
the Monte Carlo operation of the simulation program. Dropping the run from the run set . 


reduces the number of miss distaxce values that will be used to conipute statistical information 


such as mean, standard deviation and CEP; and this in turn reduces the confidence level of the ; 
statistical data. To alleviate this degradation in confidence level, additional runs are 
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*See refererc2 1 for definition of 2-cards. 
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automatically added to a run set to make up for runs terminated due to breaklock. However, 
a limit of five additional runs is built into the program to avoid ad infinitum runs (or until 
computer run time limit is reached) due to the occurrence of a very large percentage of 
breaklock flights. 


At the completion of a run set, and prior to printing out CEP data, the number of breaklock 
flights occurring, the total number of run attempts made, and the ratio of these two is 


printed out in the format shown below. 


uc MKCUCK F*BREMKLOCK FORE AKL OCK* BRE AKCOCK O49 REAKE OGK **BREAKLOCK**OREAKE DEK #* BREAKLOGKMFOREAKLOCK PB BRE Ad cOCKS- 


PURE AKLOLK® BRE AKLOCK ® 
PURCARLUCKS = —————-——. —_____ —__—— -- sa REAK LOCK 
PuREAKLOCKS THIS RUN SET HAO 17 GBREAKLOCK FLIGHTS OUT OF 32 GIVING A PROPORTION OF 5667 . *OREAKLOCKS 
© wie eae a BEAL 
cunt ARLUCK® BRE AKLOCKY 


URL ARTUCK** BRE ARCUL KO ORE AKEO COR REA RTOCK**B REAREOCK ® “BRE AKL OCK** BRE ARE OCKS* BREAKLOGK** BREAKEOCK*OBRE AKL OCK® 


4.3 Sample Run 


An example of a Monte Carlo run set and a CEP circle utilizing the OCS mode} is given 
in the following pages. The example consists of four runs out of a 25 run set that went 
into the CEP calculation. : 
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e9980RCCESCS 
FLUO CIEC3 
oS 200 CV EFGO 
eSC20GCUECOS 
sB7O0ICQESOL 
b7G ILE 
HIE IOL Cerny 
oheO0OCU ECOL 
e2OB0S OLE aL 
«26O00CQECL 
eFSCOC CEOs 
s9SC0UCCES02 
aLOOvE LEE OL 
eAe SOC OCEees 
e6L7ON0GE +O! 
oGLPOUCLEOL 
eL89IC CGE OEE 
el OFCCECL 
sSQG006GE+G2 
eSQD0GLLESO2 
e2cQ200CIE ete 
ef UOLICEsce 
oLGOCOUOCEOL 
eLQDOG CLE ele 
CTL HOG LLESOY 
eFLbOQCCE SZ 
eLIIQOVILE+Ce 
eLVGLOCVE SDS 
eULbO0gbEeas 
eBLUCOGCE G3 
eb OOO LCE Gs 
LOI IGEFGy 
LS E0IDVE+OL 
e497 QIOOECG! 
eSIOQUGCE +s 
eS5UNGCCE+COO 
665800 CUES0G 
olLSTOUGC EOL 
s7UCQICLCE+Ou 
oUF OORCL ESOL 
e7SCVOCCEG3 
el SOG QE+C2 
eLG500CCE+C2 
0835090 E-02 
ol b7 OCC E=Gt 
FOOL OIEFO2 
LF 250 0CE+0G 
eSISORCUE“L2 
eAF25UCQE+60 
eFIFCCCCEsus 
Oe 
oe 
e2H8U0UCIE+02 
AGC OUCGEFOL 
o57000CU ECR) 
oS70TUCOECCO 
o5700U0LE+OR 
oSbD NGC LEF0g 
e30L000CE+Dd 
eSH000L CEC) 
oA 00OC2TECOs 
eb2u0ulG E02 


PNA ANP Re 
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oL000090E+)1 
bl E) : 
i) 
we 3000300E0)1 
we 3O000RCE HI: 
wed OC UIICEOY 
we SQ0NI0GTHI1 
oLOOOUCOEHIL 
e LOCO eit 
we 300000 UEOI2 
Ge 
wel GQOJG0E8I4 
wel OOUIZLEeds 
wel COCIOOE+IL 
el 0GGIO0E+sL 
wel GOCIIET+31 
woh OONQONE HLS 
=Ue 
meh CL OIOCE OIL 


ar he Veawrse 


oo. bone 
0.0000 
0.0008 
=G.000 
0.0000 
0.0060 
6, 6090 
“Q.u030 
20.0090 
~0.00I9 
0.0603 
0.0060 
“0.0030 
Ge Oucu 
“Q,0ds0 
0.0050 
“020099 
“9.y979 
0.0000 
~0.4000 
=0. 0000 
“0.6080 
"0.6090 
“3.0000 
“0.0099 
=G.8050 
“0.5009 
0. 0000 
“0.6000 
=3. ouI0 
“O0.U5G8 
“G.uv00 
“Gecdiu 
“C. 0060 
ye 4060 
6.00) 
“000030 
“0.0600 
“6.0060 
0.0090 
“0. 000u 
“0-00)0 
-000930 
“0.00cu 
-0.G6000 
000090 
“Ce ¥bd0 
“06000 
CG. 6000 
Be GUlD 
3.9000 
320900 
3eu6I0 
"0.4000 
“320096 
3.4006 
360.0659 
160000 
t.u0ce 
2600)y 
1.6060 
1.C0)0 
1.6000 
70.0060 
1.q009 


* Bana 


NRE CELE PRI, AOE STN ge RE REO te TN RASTA ET 


00.6060 
“0-0.0J00 
20-0.97900 
=0"0.6090 
0-0-0000 
-0-0.000u 
=0-9.G040 
=0-9.0000 
00.6900 
=0*9.6990 
=0°0.0260 
=6-0,0u00 
-6=6.0300 
0-0. b060 
*0-0.0760 
020.0300 
-0-0.0 208 
-0-0.5900 
-0-0.9000 
-09b.0900 
-0°0.0300 
62000300 
-6-0,0000 
090.0006 
“00.000 
020.0900 
020.0300 
°0-6.0060 
20-0.0906 
“0-°0.0)36 
"0-00, 00bd 
-0-0.0009 
°0-0,300) 
°0-0,4 300 
“O-5.0900 
~"0°0.0900 
9-0, 0900 
0-0. 9000 
0-).00f0 
#00. 0029 
-0-0.0900 
“0-0,u900 
o0-0.0060 
-0-C, 0060 
*D-0.6000 
ebet. 6 ud 
-0=0.0000 
“6-6. 0000 
°0-0.0060 
oD-6.000¢ 
-0-0,0900 
“090.0000 
°0-0.0900 
-0-0,02300 
-0-9.0069 
=D- 0. 0U09 
=0-0, 0000 
=0-0.0000 
-0-6,0960 
“00.6 300 
093. 006d 
20-0.0000 
=0~0. 0000 
90-0.0y00 
0-0.0390 


ean 


CW WR W WW BWUO @BUwWWwe SG BOUWWUNaUuaecawacereocogan: 


THRST Z-OFFSET 
THRST BALPHT 
THRST OFHIE 
AUTOPILOT GY?0 
AUTOPILOT GYS9 
AUTGPILOT Grv2 
EULER ANG BPSIC 
EULER ANG ATHTO 
EULER UNG OPHIO 
JUTER UbDAL 
TNWER URAL 
TIPOFF WR MEAN 


TLAG 


TIP OFF ROLL ATE 


EZ AMPLOLAUNGH) 
FMY AMPLULAUNGH) 


EM2Z MEAN (LAUNCH) 
Fay MEAN (LAUNCH? 


NHELICUPTER VIB 


C= INOEX 
52 
» 64 
1247 
1246 
124) 
1250 
L254 
1252 
1313 
1314 
1315 
LGuL 
Laud 
1733 
1764 
1765 
1766 
1767 
$6u 
364 
362 
1754 
4753 
1752 
614 
612 
617 


eta: - weve eee ervey ee Ce Be Mwewes sts a 
43135 2 20.60 eH200GCCE“C2 -o1L 000UNICECI4 
fegi 2&2 20.60 e2SC0ICCE+G2 ef 0G0IT9EFIL 
aeC2 4 20.00 e2500C0UEtyy meL 00023 CEOIL 
$6) 2% 26.0) eSUIICCCE+CL -e20C0IOCS+)1 
$62 21 20.00 eSSUOCQUE+OL me2050,0CE+I1 
362 2 -be00 eSuGULLEECGL ect OCuaMcedt 
trou 2 20.00 eLSCOGCTEsUL ee 300CIILCEHIL 
1753 2 -0.00 ehIGIOCUESGL 23 009II T2071 
a’S2 1-660 eLvGOLCCEL =.35009C9E0I1 
d1h-G -6.60 eLS7TOCLLESGL =e FULQIIALOIL 
312-9 -0.60 eAG7OUCCE“02 ee SOUOQIDEFIL 
1738-0 -0.03 e2sCOCC&Es$2 et OCOIOC Sek 
£3230 -0.00 -.600U00CE+03 eLO0CIOVECIL 
1372-0 -0.60 eA SCOUCLEFOS o800020GE0)1 
503-0 0.U0 eA COOOCOEO1 06 4 
606-6 -0.00 2. =e 
1's3-0 +0200 e2L000CGEO2 ~e300CI00071 
1763-0 $0.60 eA GQVICEFUL =. 3000 IE+I1 
Uee=0 90.00 eAUCOQCIE+OL me JOOOIICEHI’ 
L7e5-9 0.09 eo BCUDQUCLEC2 eLGNQIVZHIL 
172-0 =d009 e2cCOGLueens e1L00LI00St31 
617-0 -O.u6 -87UCULUE=83 eLUOG000S+I1 
d£3-G 0,60 eB7UI0CSE-03 eLO0QIIIEMIL 
917-0 -0.60 eSLOICCCE=04 =o 3O9fVIGE+IL 
o23-C -0.60 eSCOOOCEE08 =e 3UQUUOCEFIL 
312-0 90.60 oS5°50GLS E04 el OOCIIUTHIL 
d1L°0 -0.08 eo 5uGOCCKE=61 e1000I0IEFI4 
4315-0 -G.00 -.1500CCLE+U5 eLOVGINOE+T1 
755-6 30.00 =. 20C0.6C E01 oLOCCSI9E4I1 
AP E3+9 90.00 meds TO0CIE+CL 010603) 7E4)1 
1315-0 -0.00 -.G20000CE+04 eLQQ0INGE+) 1 
1753-0 -0.60 -eb9CCOCLEFCE eiCO0ICGEHOL 
1515-20 $0600 =. 7K000CGE+O4 e1 D00I00ECI1 
3235-6 -O0ad el. COOCVE+CA -de 
“0-0 -0.50 7. we 
MONTE CARLO INTTIAL CONDITIONS 
MC-VALUE MEAN DISTRIBUTION LOWE 
8.2268072 0,000)000 NORMAL 
236.561 0.000 UNIFORM 
me18) C000 JNIFORM 
2378 0.006 UNIFORM 
=2329 O.duu ¢ UNIFORM 
72174 Go0C9 UNIFORM 
-0377 GeOul JEIFORM 
0c8h 0,000 UNIT FORIA 
we tue Good UKIFORM 
~.0us Q,000 UNIFORM 
--CC3 v.00 VNT FORM 
0 C74 0.065 UNIFORM 
=2220 Ge0LG UNIFORM 
8eL197569 2640002909 NORMAL 
260645006 Q.0G65902 NORMAL 
-.069%272 Q.U003303 NORMAL 
ou 605035 O.80G)95] NORMAL 
00445268 0.6000009 NCRMAL 
1.775 0.060 UNI FORM 
2597 0.0Cd UNZ FORM 
22924 0.000 JNIFORM 
«858 C.006 Uni FORM 
“4.706 4.960 UNIFORM 
4.389 0.006 JAI FORM 
0 885407 eOSC2)Uy NORMAL 
ou 379133 o60S5705 NORMAL 
06008617 00008703 NCRHAL 
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aveves 
4.0030 
120090 
12G6C)0 
1. GO00 
4.06)0 
4.0090 
3. 0000 
3.6070 
3. 0600 
3.6030 
3.6009 
-C. Gu000 
“0.0000 
9.6000 
=0,00)0 
6.0669 
3.00Cu 
3.6900 
3.0030 
“0.0000 
=0,0099 
“0- a0ua 
=0.Guv0 
3.0000 
3. GU) 
on. 2906 
-G.cCtC 
-0.00)0 
0.0050 
Ce 00)) 
“C. aud 
0.0020 
“0.0000 
70.0008 
=G. 0000 


R BOUND 
~-3.000 
0.090 
=1.000 
71.000 
=1.600 
21.600 
-1,000 
-1.090 
~1.000 
°1.090 
wiL..00 
-1.690 
1.000 
~3,.600 
-3.090 
~3.000 
-3.000 
-3.090 . 
71.090 
“1.000 
74.000 
= 3.000 
~3.090 
=3.050 
-3.UN0 
~3.090 
=3.000 


-C20.G090 
-0-0. 9000 
-9-0.00.0 
=6-0.0000 
-o-C.COC0 
-0°0.0369 
-O-6.0006 
-0-2.6000 
"0-0.0000 
“070.0060 
=0-0,05C0 
-0-0,0000 
60.0300 
00.0000 
-0-0.0000 
-0-6.0000 
*0°0.0u00 
“Ge-(.0J300 
*0-0.0080 
-G-0. 03550 
-0-0.00C0 
“0-0.0000 
-b-0.0GU0 
“u-0.5I00 
06.0680 
“020.0192 
-0-¢.5000 
-C-0.03C0 
*0-C.09Cd 
-0-6.0000 
-0-6.0000 
-0-0,0)59 
*0-0.00u0 
-0-0,0000 
-0-0,.0000 


UPPER BOUND 
3.000 
350.0060 
1,00) 
1.00) 
1.00 
1.050 
1.000 
1.00) 
4.069 
42003 
1.053 
4.09) 
2.000 
3.02 
3.0u9 
3.0u9 
3.609 
3.099 
1.00) 
41.04) 
1.004 
3.099 
3.009 
3.099 
3,209 
3.04] 
3.009 


ie ee ee pe en ma 2 


seep tes one 


Tinea oN92; 0G STE? Sizes 2.0006Cu0E-03 aie 
FRONT LUG SLZARS ANZ. T2 a.60E-92 REL VEL # 2sL18¢C2 PLIC4 MOMENT @ -5.45266L 
RANGO = 3.5868 


camaro Meet cee 


TIPOFF WATES=-ROLL 2 SL PITCH 2 “4.3 YAW ® o3el 
REAR LUG CLEARS ZATL T 2 eAL2CREL VEL = 232.928 RAIL FORCE = 06063 


RANGO = 569739 
TIMER o112290G SEP SIZE 8.350 UL00E-03 


GURNOUT TIMEa 3.0041 3ES. 


EL EE EE I 


TERY UIP USI PererVlie ir evivinirririr iii rrirvirirrrrirerrisvrrrrivriniriirtliviriiirirririirryig) 
OCS BLIND RANGE SIGNAL HOUU AT TIME = 6224 RANGE 2 329279 


BREE SHUTESRFREST GHSEHFFSF RE SSHE ES HAE RES SEE SEHSEREFESSHSSESESEL ESE PSST EFF RE FEVER ECR SEEAAER OPER SERRE 


FTE TT Rae TY 


o++MAX BREACLOGK VALUE = 0283)? IN. PITCH 


Ob+HMAK BREACLOCS VALUE 027224 IN YAW 
RUN NUMBER = 1 


UISS OLSTawCE = 4.752491 25E400 


FLIGHT TIME = 6053603 95E+GU : 
ROELX © 6.73 31478E-01 ROELY = &,514565755*31 WELZ 2 4 7 9539GE NO 
. RYFP mm .-5H5K947E5)4 QzFP =e 72962098499 
82 Bese Sy Sages es 
RY FERS SE Sassi gs 
i 
i 
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a 
1 


ca ram s2 soil daciiecia i Saal 
“ iis Naa ch al Pas ee te ee eS SEE DRU TIENT Y SORTER ETRE WOVE U Te Fone 
coed hab stage ee afte ne Re eet ee att la men Bn aed haa asin malate la 


TaPJT OATAa 


& -0-6 - 


C-IN 

52 

Si 

1247 

1248 
1249 

1250 

1251 

1252 

15i3 

1314 

1315 

1461 

1462 

41733 

1764 

1765 

1766 

1767 

365 

361 

=a. eget, 362 
1754 
1753 
1752 
611i 
612 
617 
618 
sree siz 


TIMia 0027000 


FRONT LUG CLEARS RAIL 


REAR LUG CLEAIS ATL 


TIMES e112.900 SfeP Siz 


QURNOUT TIMZS 340021 3&5. 


COSHH SEP ERT SS SETH CHSEHEE SSE HFETSSSSEEHS ESOS SSSA GSES SE SSSHR FEARS EUSSHS PERE SIL SECIS GSES FHAPSR ESTEE EREEARED 


OGS BLIND RANGE SIuhAc HOLD A 


SS ee ee ee ee ee 


-0-0.0900 
LOWER BOUND UPPZR BOUND 

3.006 3,090 

G.ud0 350.0u0 
-1.606 1. bud 
-1.00C 1.099 
o1.00¢ 16099 
“1.3390 1,02 
-1.C600 12040 
-1.096 1.000 
-1.00C 1.0659 
“1.000 1.009 
=1.000 1.00) 
71.099 1.009 
-1.0C6 1.090 
3.000 3,009 
-3.000 30 Oud 
=3.000 3.090 
-3.65C 3.000 
=3.900 32930 
“1.600 te Iu 
-1.690 1.003 
“1.000 1.069 
-3.000 3.053 
-3.006 3,003 
3.000 3.069 
-3.000 3.099 
=3.604 3.900 
=3.u00 3.009 
“3.006 3.0y0 


. 


PKTC4 MOMENT 2 -5.43E4CL 


36.3 PITCH = 


G.u0 =C. “Ge “0.0030 
MONTE CARLO CNITIAL CONDITIONS 
DEX MC- VALUE MEAN IISTRIDUTION 
269129063 0.0003009 NORMAL 
163.976 0.000 UNTFORN 
me K7& 6.000 UNIFORM 
=oCc0 G.00¢c UNTFORM 
Car | €.090 UNS FORM 
ee 315 0.000 UNIFORM 
0180 6.006 UNZ FORM 
70223 G.g0c JNIFORA 
23 0.06C - UNIFORM 
wel Ok 0.0u0 JNIFORM 
052 C.0CG UNIFORM 
22152 0.009 UNIFORM 
abl 0.000 UNIFORM 
36.2631223 20.0060)90 NORMAL 
eutesurs O.0090090 NORMAL 
01738540 0.000239) NORMAL 
~2 eb 695206 0.0009909 NORMAL 
=e lb 13949 €.0009002 NCRMAL 
720606 0.000 UNIT FOR 
1.696 0.000 UNIFORM 
6713 C.060 UNIFORM 
0389 0.000 UNIFORM 
=3.681 =4.90C UNLFORM 
1.190 0.000 UNI FORM 
0U418523 06502990 NORMAL 
09615168 20503000 NORMAL 
e0LC3060 eOG0370yG = NGR MAL 
0 UCC8456 00008762 NORMAL 
Ez 2.606 du QUE-03 
T = 865E=62 REL VEL @ 2eL1E*b2 
RANGO = 3.5829 
TIPOFF RATES“-ROLL 2 
Tes eLL2UREL VEL = 232.875 RATE FORCE 
RANGO = 5.9769 
E= 8.350 5000E-03 


T TIME = 6.624 RANGE = 305.73 


~41.32 


“21.3 YAd =. 


SEOPER SSH HS ETRE SESETHS SHH STHETSESESE SSP TE RFS SVEFE EEE EES SETAE EESES EKER SESE SEER EIS SERPTSPHESERE RE RES 


RUN NUNDER = 2 


*++MAX BREAKLOCK VALUE = 


*++MAX BREACLICK VALUE = 


MISS DISTANCE = 368s 27765E+00 
FLIGHT TIME = 6 546726 7E900 
ROELX = -Se72U94BLE*O4 RDELY = 5e4297058E-31 
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RDELZ 


0389)9 


=234518 


S-7595870E+OG 


71.9 


InN PITCH 


In ¥ 


An 


{ 


Ess’ r oe EY Se NE OR REE EP PT ELIE TS I TE 


BELT FLL EGRET TL EEE TN NS NT ELT ee eT 


r 

Be 
ba 
\ 


RYFP Se 32KSLLOE3L RwFP ok 3.8057098E+00 
INPJT DATA 
Be eee 20-0 -0.00 <0. --- 1+ Qe “0.0000 =-0-0.0380 , “ 
HONTE CARLO INITIAL CONDITIONS 
Ce INDEX MCG-VALUE MEAN JISTRIBUTION LOWER BOUND UPPER BOUND 
52 5001629962 02000320) NORMAL -3,000 32009 
: 54 248.545 0.000 INTFORM 06260 350.009 
some a+ L267 : e193 Ge000 UNIT FORM toe =1.000 1.009 
1246 «066 0.000 UI FORM 4.600 22009 
1249 -0118 6.000 JINIFORM =41.000 4200) - 
1250 2261 0.000 UtiT FORM -1.000 1.069 
4251 22f89 0.090 UNIFORM “4,000 2-90) 
41252 2163 0.000 UNIFORN =1.600 1.000 
tee tee eee ee ~ £343 wo ws + C03 6.000 UNIFORN -1,000 1.009 
1314 otC3 G.600 JNIFORM ~1.000 1.609 
£315 , ef G3 C.000 UNI FCRM -1.690 1.00) 
L464 0305 0.000 UNIFORM 71-000 4.000 
Sees 14020 : - 0156 0.060 JNE FORM ~1.000 1.093 
1733 G3.34615734 26.0009900 NORMAL 3.000 3.008 
a . 7 te - + 2766 tote 09046708 Oe0G0I702 NCKMAL : -3.00 3.002 
4765 -e6206029 0.600030) NORMAL -3.600 36000 
1766 ~ 2b 677724 060900063 NORMAL -3.006 - 3-00) 
1767 oF 243761 620000063 NORMAL -3,000 3,00) 
360 2538 0.000 UNIFORM =1,000 1.000 
: 36. 2.257 0.000 URI FORM -1.000 1.000 
RENEE oo Hes ards aermins «patie 7 - + 362 Secs 0675 0.000 URIFORK - -1,.600 4.660 
1754 0855 0.960 UNIFORM “3.000 32 Ou0 
1753 o3.641 -4. 900 UNIFGRM =3.600 Je GLO 
1752 -2.797 0.000 JSNIFORM 3.600 3. Gul 
- 614 : 9665639 005000G0 NORMAL 3,000 3. 000 
612 6585906 00503067 NORMAL -3.00C 3. Gud 
617 “ 09663758 00053792 NORMAL -3.000 32°20 
613 00008585 0006870) NORMAL 3.660 , 3.090 
TIMES - 0602. 5906 SiEP SIZE2 2eCOOCLQUE-03 
FRONT LUG CLEARS RAIL T= 8.666-02 REL VEL = 2eL1E+02 PITS4 MOMENT = -6.4)E¢UL 


RANGO = 3.5828 
wy hale: vies Tees a3 ee TIPOFF RATES=-ROLL &- 4304 PITCH = -22.) YRA = bee 


RANG 3 5.976% 
TIMES . 04125900 SIEP SIZE= 6-350 TUORE*02 


REAR LUG CLEARS RAIL Ts eAL2GREL VEL = 232.375 RAIL FORCE = “43.56 


BURNOUT TIMES J.0014L 5E5. 
SORE EE ERED ESSERE SEEE EO EDO REDD SOORUSESES TR ERSEER ERED EERESSE SECU ESE RUS HEE AE RE EERE OUURE EEE EY EERE MEEE ES 


OCS BLINO RANGE SIGNAL HILO AT TIME = 6.16 RANGE 2 305494 
COON OF TTR SR EEEE SEES FS BEE EERE OE ERUERE AEE EREEEL SUSUR EEEPEFEEEU ERE REE SER OSE ERE RRR ERIE EERE HERE BUSES 


+OOHAX BREAKILOCK VALUE = 436498 IN PITCH 
CeOMAX BREACLOCK VALUE = 422691 IN YAW 
RUN NUMER #3 
MISS UISTANCE =  5,9798868E+00 


FLESHT TIME = 6,4) 6¥555€+00 
RDELX © *6,5260502E=01 ROELY = 5.2700933E-214 RDELZ = 389041708490 
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cpedstn dee tera aspirin Lee ducae seen Scalise siatdalil 


sedate nt nbd ian SLS Ek AS a nasi na cai Nh tad is tina ited aint edn Nemabbldanina sor 
cea kan re te eretat ea 


eke re awa ae pated oe aE ae rie 


2 Se EE Em res ce eo pot ee eater eee cme age aga er 


RYFP GeD70357KE-)E RZFP 2 3.9487081E+00 


INPJT OATA 
Gos : 23-0 -0.60 <9. nde =0.0030 -0-0.0000 
MONTE CARLD INITIAL CONJITIONS 
C-INOEX MG VALUE MEAN JIS TRIBUTION .OWER BOUND UPPER BOUND 
52 ~30.7594857 0.009022) NORMAL =3.000 3.009 
St 20.578 6.006 UNIFORM 0.090 360.093 
12407 2352 92000 UNI FORA 1.600 10060 
1248 2610 4.006 UNIFORM 1,600 1.000 
1249 2195 4. 0C6 JNIFORM =1.066 1.903 
1250 2223 6.060 UNIFORM 1.000 1-003 
Byide B- fy eR oe ote 2251 -.073 0.000 UNIFORM 1.606 1.0ua 
‘ 1252 2150 0.000 YNLFORM 1,090 1.009 
= sd - 1313 -.bv2 0.009 UNIFORM -1,.606 1.099 
L314 ofOL 0.2060 JNIFORK *=1.000 1.Cu3 
: 2315 3663 0.000 UNTFORN -1.006 4.669 
. Loul ole C.600 JNIFORM -1.006 212009 
: df 72 142 2.043 0.0ce0 JShIFORM 1.000 1.000 
1738 204428024 20.0600362 NORMAL -3.Cu0 3.03) 
em = L764 01296747 =§.0000)00 NORMAL -3.000 3.03) 
1765 06250585 0.0063000 NORMAL 3.000 3.0699 
1766 ; 20722963 G.00002700 NORMAL 3,600 3.000 
: is 1767 - 20855807 0.0603000 NORMAL -3.000 3.900 
: 366 o3il 4.000 UNLFORM #1.600 1.099 
; Jor 2.630 6.000 UNIFORM 1.000 4.009 
: ve hw aed 362 2.684 6.006 UNIFORM 1,090 1.000 
: 1754 + 366 0.600 UNIFORM 3.006 3.009 
$ . - £753 2.657 4.960 UNIFORM -3.600 3.000 
‘ i752 20346 0.000 ’ UNIFORM = 3.000 3.909 
: . 611 0v693596 2G5C39CG NORMAL -3.000 3.099 
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i TIPOFF ZASES=-KULL = 4205 PITCH = =i60d YAS # 1.55 - 
i “ REAR LUG CLEAAS RAIL T= eALZUREL VEL = 2326525 RAIL FORCE ®© -35,63 
\ RANGO = 5.9740 
TIME® - 24125900 SFEP SZZEs Be I5CCUOGE-03 


BURYOUT TIME= 3.0612 5ECe 
POVYOU OTOP TPT ie Tri Tiere er irrerirvrr ete rrerrrevirvrivisrirriservirevrarertrirtrvrririrerrirrivirr! 


OCS BLIND RANG? SIGNAL HILO AT TIME = 6244 RANGE = 316.34 
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+eeMAX BREAKLOCK VALUE = 35553 IN PITCH 
+e+HAX BREAKLUCK VALUE = 027738 IN YAW 
RUN HUMBER = & 
“MISS DISTAKUE a wo 4) 723046400 
FLLGHT TIME 2 6.775467 3E006 j 
RDELX = -5,1798049E-04 RUDELY 2 -Gel4H4651E-34 QELZ = HL 3SA]VOTMGEFOC { 
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THE 99 PER CENT CONFIDENCE CIRCLE RADIUS IS 14418 
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5.0 COMMENTS 
5 
| 5.1 Integration Synchronization With Sample Period 


t 
Numerical integration must be synchronized with sample period (r @ ZOH) in order to insure 


accurate integration. Logic is built into the seeker subroutines to insure that integration and 


Ee ea Re ee ey ee a ee et 


sample period are synchronized, This is accomplished by: 


At = T/AINT(T/At, orp)! 


where 


T ' « sample period (SEC). 


4 
: 


4tunput — input integration stepsize GEC). 
AINT (X) — computer center library function that integerizes 
the argument (X). 
At — computer integration stepsize that the program will use. 


The above function will always compute an integration stepsize that is equal to or greater than 
the input stepsize, Since there is an upper bound on the stepsize that can be used to integrate 
the differential equations in this simulation program, there is the possibility that a ‘stepsize 
larger than the upper bound will be computed. (Upper bound is approximately 12.5 millisec, 
with the exception of the OCS seeker model S2 which has an upper bound of approximately .5 
millisec), Therefore, one should insure that a reasonable stepsize is input and verify that a 
reasonable stepsize is computed. For example, if the sample period is 16.7 millisecs, then 
a stepsize of 8.35 millisecs or less must be input to insure that the computed stepsize is 12.5 


millisecs or iess. 


* A discussion of errors caused by numerical integration of sample data is given in 


Reference 5. 
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APPENDIX 


Monte Carlo 6-DOF Program Listing 


A FORTRAN IV listing of the Monte Carlo 6-DOF Program that is operational on the 
MICOM CDC 6600 computer, SCOPE 3.4 operating system, is given in the following 


pages. 
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